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Abstract 

We consider a modified gravity fluid on a Randall-Sundrum II 
brane situated at y = 0, the action containing a power a of the scalar 
curvature. As is known from 4D spatially flat modified gravity, the 
qh| presence of a bulk viscosity may drive the cosmic fluid into the phan- 

tom region {w < —1) and thereafter inevitably into the Big Rip singu- 
bl), larity, even it is initially nonviscous and lies in the quintessence region 

{w > —1). The condition for this to occur is that the bulk viscosity 
contains the power (2a — 1) of the scalar expansion. We combine this 
with the 5D RS II model, and find that the Big Rip, occurring for 
Q^ I a > 1/2, carries over to the metric for the bulk metric, \y\ > 0. Ac- 

t^^ ' tually, the scale factors on the brane and in the bulk become simply 

proportional to each other. 
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r^ • 1 Introduction 

H 

5r ! Modified gravity theories in 4D continue to attract interest; this obviously 

being related to observations, for instance the measured redshifts from type 
la supernovae [H [2[ [3]. The data may be reconciled with the concept of dark 
energy, with a cosmic fluid with a complicated equation of state, or with a 
scalar fleld having quintessence or phantom behavior. An extensive recent 
review is given by Copeland et al. [4J . The equation of state for the cosmic 
fluid is conventionally written as p = wp, where w = —1 corresponds to a vac- 
uum fluid (cosmological constant), — 1 < w < —1/3 to a quintessence fluid. 
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and W7 < — 1 to a phantom fluid, having the bizarre property of predicting a 
Big Rip singularity in the future. 

In the present paper we combine essentially two kinds of theories: 

1) First, we assume the modified 4D action in the form given by Eq. flTSl) 
below. The integrand contains a power a in the scalar curvature. This 
model has recently been studied by Abdalla et al. [5J. As a generalization, 
we include a bulk viscosity in the fluid being proportional to the {2a — l)'th 
power of the scalar expansion. Cf. Eq. (13T|) below. This viscous model has 
been studied repeatedly in the recent past [6l [71 [H [9] . A striking property of 
this kind of theory is that it leads to the Big Rip singularity, if a > 1/2. 

2) The next step in our analysis is to combine the above 4D theory with 
the 5D Randall-Sundrum II model, [TD], where there is a spatially flat brane 
situated at |/ = 0, surrounded by an AdS space. Our main result is to show 
how the mentioned Big Rip singularity on the brane becomes transferred to 
the bulk: the scale factor away from the brane {\y\ > 0) becomes actually 
proportional to the scale factor on the brane itself. In this sense the 4D and 
5D gravity theories are closely intertwined. In this respect there is no essen- 
tial difference between the modified gravity and the Einstein gravity. They 
behave qualitatively in the same way, only with the characteristic difference 
that the strength of the singularity becomes larger for increasing values of a. 

2 Basic formalism 

Assume, as mentioned, that there is a spatially flat {k = 0) brane located 
at the fifth dimension y = 0, surrounded by an Anti-de Sitter (AdS) space. 
If the five-dimensional cosmological constant A(< 0) is different from zero, 
this model is the Randall-Sundrum II model (RSII) |lDj. We shall take the 
metric to have the form 

ds"^ = -n'^{t, y)dt^ + a^{t, y)6ijdx'dx^ + dy'^, (1) 

where n{t, y) and a{t, y) are determined from Einstein's equations 

Rab - ^gAsR + QabJ^ = i^^Tab- (2) 

The coordinate indices are numbered as x"^ = (t,x^,x^,x^,y), with k^ = 
SttGs the five-dimensional gravitational coupling. Einstein's equations in a 



coordinate basis with the metric ([T]) have been given before (cf., for instance, 
Refs. [m 1121 [ISl El [ISl [IS] , but for convenience we give them also here: 
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Overdots and primes mean derivatives with respect to t and y, respectively. 
As the bD space outside the brane is taken to be empty, the components of 
Tab are different from zero only on the brane. 

Consider next the form of Tab- Let U^ = {U^, W) (Greek indices /U, z/ e 
[0,3]) be the fluid's four- velocity on the brane, and let a denote the brane 
tension, assumed constant. Moreover, let h^u = g^v + U^Uu be the projection 
tensor, and p = p — 3Ho( the effective pressure with Hq = do/ao the Hubble 
parameter on the brane and ( the bulk viscosity. The shear viscosity is 
omitted due to the assumed spatial isotropy. 

It might be noted here that a viscous fluid may be considered as a specific 
example of an inhomogeneous equation-of-state fluid introduced in Refs. [171 
HE]. 

As gauge condition we take no(t) = 1, this meaning that the proper time 
on the brane is the same as the cosmological time. The energy-momentum 
tensor can accordingly be written as 



AB 



;i, 0,0,0), and let a 



We shall work in the orthonormal frame where If^ 
subscript zero refer to the brane. 

The junction conditions aXy = have now to be taken into account. They 
express that the metric is continuous across the brane, but its derivatives are 
not. From Eqs. ([3]) and (jl]) we get, for the distributional parts. 
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[n'] = -K\-a + 2p + ?,p), 
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where [a'] = a'(O^) — a'(0 ), and similarly for [n']. For the nondistributional 
parts we get 
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Assuming no energy flux to occur from the brane, we have Tty = 0. It implies 
that 

n{t,y) = ——r- 12) 

ao{t) 

for arbitrary y. Then from Eq. (ITOl) we get, upon integration with respect to 

y, 
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Here C = C{t) is an integration constant with respect to y. The C term 
is called the "radiation term"; as it is not of main interest here, it will be 
omitted in the following. 

Now setting y = we get for Hq = do/ao 
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Recall that A and a are constants, while p = p{t). The generalized Friedmann 
equation (iMl) can be contrasted with the conventional Friedmann equation 
in four-dimensional space (still with k = 0), 
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with k| = 87rG4. The essential new feature of Eq. f lT^ is thus the occurrence 
of a p^ term. 

Let us observe the solution for a^it) from Eq. flT^ if p = 0: 
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where 

A=1a4.V, (17) 

Co being a new integration constant. The scale factor is thus exponentially 
increasing, qualitatively as in the ordinary de Sitter case. 



3 Modified gravity on the brane 

In this section we consider the fluid - Einstein or modified fluid - on the 
brane y = 0. We shall derive how the Hubble parameter H varies with time 
t, leading eventually to the Big Rip. 

We adopt the following 4D gravity model: 

S = ^ I d^x^TdU^R^ + L^), (18) 

where /o and a are constants [a may in principle be negative), and L^ is the 
matter Lagrangian. This model has been studied before; cf., for instance, 
Refs. [5|, [71 [8]. The case /o = 1, a = 1 corresponds to Einstein's gravity. 
(More complicated /(-R) theories have been discussed at various places, for 
instance in Refs. [T71 [HI [19] . We may also mention here that a general review 
of modified gravity can be found in Ref. [20], and recent reviews of viable 
/(i?) gravity theories unifying dark energy, inflation, and dark matter, can 
be found in Refs. [HI [22].) 

The equations of motion following from the action (TT8|) are 

— -;^fo9tJ.uR°' + afoRfiuR"' 

- a/oV^V,i?"^i + afo ^^.V^i?"-^ = ^T^,, (19) 

where T^jy corresponds to the term Lm in the Lagrangian. (The proposal of 
taking f{R) on the brane was considered also in Ref. ^23j.) 

We take the equation of state for the fluid to have the conventional form 

p = wp={j- l)p (20) 

(more complicated forms for the equation of state have recently been inves- 
tigated by [T71 nil [21]). li w = —1 OT p = —p we have a "vacuum fluid", 
with bizarre thermodynamic properties such as possibly negative entropies 



As is known, cosmological observations indicate that the present uni- 
verse is accelerating. Moreover, based upon the observed data it has been 
conjectured that w is a varying function of time. For instance, as discussed 
in Ref. |26], w might have been around at redshift z ~ 1 and may be 
shghtly less than -1 today. Perhaps is w even an oscillating function in time. 
Recent discussions on possible forms of the equation of state are given, for 
instance, in Refs. [27| and [2B]. In view of these circumstances the analysis 
of a possible crossing of the phantom barrier w = —1, from the quintessence 
region (—1 < w < —1/3) into the phantom region w < —1, is obviously of 
physical interest. It ought to be noted that both quintessence and phantom 
fluids lead to the inequality p + 3p < 0, thus breaking the strong energy 
condition. The Big Rip singularity has been discussed in various papers; cf., 
for instance, Refs. [29], |30l [3T[ [321 [33] . Specifically, the future singularities in 
modified gravity were considered also in Ref. [34] . 

Of main interest is the (OO)-component of Eq. f[T9l) . Observing that i?oo = 
—3a/ a, R= 6{H + 2H'^), as well as Tqq = p, we obtain 

^/oi?" - 3afo{H + H'^)R'-^ + 3a{a - l)foHR"-'^R = kI p. (21) 

An important property of Eq. f[2Tl) is that the covariant divergence of the 
LHS is equal to zero [35] . 

V%. = 0, (22) 

just as in Einstein's gravity. Energy-momentum conservation is a conse- 
quence of the field equations. This leads to the energy conservation equation 

p + {p + p)3H = %H\ (23) 

We now differentiate the expression f[2Tl) with respect to t, and insert p from 
Eq. f[23|) . After some calculation we obtain 

^7/o^" + 3afo[2H - 3j{H + H^)]R^-^ 

+ 3a{a - l)/o[(37 - l)HR + R]R''-'^ + 3a{a - l)(a - 2)foR^R''-^ = 9kICH. 

(24) 
Recalling that R = 6{H + 2H'^), we see that this equation is a complicated 
nonlinear differential equation for H{t). It is best discussed in terms of 



examples. We are interested in solutions that are related to the Big Rip. We 
shall look for solutions having the form 

TT 

H=^, where X = l- BH, t. (25) 

Here H^, is the Hubble parameter at present time t = (usually called Hq 
but we are reserving the subscript zero mainly for brane entities), and B is 
a nondimensional constant. For Big Rip to occur, B has to be positive. 

3.1 Einstein's gravity fluid 

As mentioned above, this case corresponds to /o = 1, a = 1. As for the bulk 
viscosity, we shall take ( to be proportional to the scalar expansion 9 = 3H 
through a proportionality constant, here called Te, 

( = te9 = SteH. (26) 

This form is of particular physical interest. Namely, as shown in Ref. [6], if 
Te is large enough to satisfy the condition 

X = -j + 3kIte>0, (27) 

then the equations of motion lead to the Big Rip singularity in a finite time 
t. Even if one starts with a state where the fluid is nonviscous and lies in 
the quintessence region (7 > 0), the imposition of a sufficiently large bulk 
viscosity will drive it into the phantom region and thereafter inevitably into 
the Big Rip. 

From the governing equations we now get 

B = Ix, (28) 



H, = \jl4p*, (29) 

Pe = j^, (30) 



where p* is the t = value of the energy density. 



3.2 Modified gravity fiuid 

Assume now that /o and a are arbitrary. Let the bulk viscosity for the 
modified fiuid be denoted by Ca- As in Refs. [H [9] we model Ca by setting it 
proportional to the (2a — l)'th power of the scalar expansion: 

C = rj''^^^ = r„(3iJ)2"-i. (31) 

The main reason for this assumption is that it fits well with the governing 
equation for H as well as with our previous assumption (12^ : the time- 
dependent factors in Eq. (^^ automatically drop out, and we remain with 
the following equation determining B: 

(B + 2)""^{9(2 - a)7 + 3[a + 37 + a{2a - 3)(37 - 1)]B 

+ 6a{a-l){2a-l)B'} = ^ Q\^. (32) 

This equation is complicated, and is best discussed in terms of examples. For 
instance, if a = 2 and 7 = (the latter condition corresponding to a vacuum 
fiuid), then Eq. ([32D yields the following cubic equation {to — >■ T2): 

B' + 2B' = ^. (33) 

0/0 

There exists one single positive root of this equation, as long as the RHS is 
positive. This root is caused by viscosity, and leads to the Big Rip. 

We note also the general equation for B following directly from the energy 
conservation equation fl23|) for the modified fiuid. 



Pa + (Pa + Pa)^H = 9CaH^ (34) 



namely 



where we used 



37 3r,(3^)^ 

(a = Ta [-^J . Pa = ^- (36) 

For simplicity we have assumed the same initial conditions at t = for the 
modified fiuid as for the Einstein fiuid, viz. p*Q, = p^E = P*, H^:a = H^.e = 

8 



4 Implications for the 5D theory 

We are now equipped with the necessary background to see how the modified 
fiuid on the brane effects the 5D brane physics. Consider first Eq. (fT^ on 
the brane (recall that this is a 5D, not a 4D, equation). It is natural from a 
physical point of view to use the expressions for p{t) from the previous section 
as input quantities in this equation. Comparison between Eqs. f l36|) and fl30l) 
shows that we can regard pa = p*/X^" as a generic equation common for 
the two cases, only with a = 1 in the Einstein case. For the 5D scale factor 
ao(t) on the brane we obtain thus 
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As we shall be mainly interested in the behavior near Big Rip, we consider 
times close to the singularity time tg = 1/{BH^), where we get approxima- 
tively 






p* 
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The quantities A and a, characteristic for 5D theory, are here neglected. The 
solution of this equation is of the form 

(kV6)p* 



ao(t) ~ exp 



(39) 



_(2a-l)(5/7,)2-(t, -t)2"-i_ 

The scale factor on the brane has thus an essential singularity at t = t^, if 
a > 1/2. This behavior incorporates both the Einstein gravity, and the R^- 
modified gravity (a = 2), considered in the previous section. The singularity 
is stronger the larger is the value of a. Moreover, the divergence is stronger 
than the power divergences found for viscous 4D cosmology with account of 
quantum effects |36]. If a < 1/2, oq does not diverge at tg. Note that there 
is a relationship between p^ and if^, in the expression (l39l) as following from 
Eq. ([H]) taken at t = 0, 
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Now return to the bulk case, considering Eq. flT3l) for arbitrary y. 
C = as assumed (recall that also /c = 0), we obtain as AdS solution 



(40) 
When 
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(41) 



with /i = J— A/6. And this brings us to the following important conclusion: 
The Big Rip divergence on the brane, present as we have seen when a > 1/2, 
becomes transferred into the bulk. The bulk scale factor a{t, y) diverges for 
arbitrary y at t = tg, ii ao(^) diverges at tg- This result could hardly have 
been seen beforehand, without calculation. There is moreover no particular 
difference between an Einstein fluid and a modified gravity fluid in this re- 
spect; they behave essentially in the same way. It may also be of interest to 
note that if the brane is tensionless, o" = 0, then the bulk solution becomes 
quite simple, 

a\t,y) = ^al{t)[l + cosh(2^y)]. (42) 

The bulk scale factor thus increases exponentially on both sides of the brane, 
for large \y\. 

It would be of physical interest to understand the simultaneous occurrence 
of singularities on the brane and in the bulk - perhaps there are quantum 
effects at play here. 



10 



References 

[1] A. G. Reiss et al, Astron. J. 116, 1009 (1998). 

[2] S. Perlmutter et al, Astrophys. J. 517, 565 (1999). 

[3] J. L. Tony et al, Astrophys. J. 544, 1 (2003). 

[4] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 
1753 (2006). 

[5] M. C B. Abdalla, S. Nojiri and S. D. Odintsov, Class. Quant. Grav. 22, 
L35 (2005). 

[6] I. Brevik and O. Gorbunova, Gen. Rel. Grav. 37, 2039 (2005). 

[7] I. Brevik, O. Gorbunova and Y. A. Shaido, Int. J. Mod. Phys. D 14, 
1899 (2005). 

[8] I. Brevik, Int. J. Mod. Phys. D 15, 767 (2006). 

[9] I. Brevik, Gen. Rel. Grav. 38, 1317 (2006). 

[10] L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999). 

[11] I. Brevik and A. Ballanger, Phys. Rev. D 69, 024009 (2004). 

[12] P. Binetruy, C. Deffayet, U. EUwanger and D. Langlois, Phys. Lett. B 
477, 285 (2000). 

[13] P. Binetruy, C. Deffayet and D. Langlois, Nucl. Phys. B 565, 269 (2000). 

[14] I. Brevik, K. Ghoroku, S. D. Odintsov and M. Yahiro, Phys. Rev. D 66, 
064016 (2002). 

[15] I. Brevik, K. B0rkje and J. P. Morten, Gen. Rel. Grav. 36, 2021 (2004). 

[16] I. Brevik, J.-M. B0rven and S. Ng, Gen. Rel. Grav. 38, 907 (2006). 

[17] S. Nojiri and S. D. Odintsov, Phys. Rev. D 72, 023003 (2005). 

[18] S. Capozziello, V. Cardone, E. Elizalde, S. Nojiri and S. D. Odintsov, 
Phys. Rev. D 73, 043512 (2006). 



11 



[19] K. A. Bronnikov and S. G. Rubin, Gravitation & Cosmology 13, 1 
(2007). 

[20] S. Nojir i and S. D. Odi ntsov, Int. J. Geom. Meth. Mod. Phys. 4, 115 
(2007) [hep-th/06012l3]. 



[21] S. Nojiri and S. D. Odintsov. [aFXiv:0801.4843l [astro-ph]. 

[22] S. Nojiri and S. D. Odintsov. ' i?Xiv:0807.0685i [hep-th]. 

[23] S. Nojiri and S. D. Odintsov, Gen. Rel. Grav. 37, 1419 (2005). 

[24] R. A. Sussman. [i[FXiv:0801.332 4 ] [gr-qc]. 

[25] I. Brevik, S. Nojiri, S. D. Odintsov and L. Vanzo, Phys. Rev. D 70, 
043520 (2004). 

[26] A. Vikman, Phys. Rev. D 71, 023515 (2005). 

[27] I. Brevik, O. G. Gorbunova and A. V. Timoshkin, Eur. Phys. J. C 51, 
179 (2007). 

[28] I. Brevik, E. Ehzalde, O. Gorbunova and A. V. Timoshkin, Eur. Phys. 
J. C 52, 223 (2007). 

[29] R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 
91, 071301 (2003). 

[30] B. Mclnnes, J. High Energy Phys. 0208, 029 (2002). 

[31] J. D. Barrow, Class. Quant. Grav. 21, L79 (2004). 

[32] S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71, 063004 
(2005). 

[33] S. Nojiri and S. D. Odintsov, Phys. Lett. B595, 1 (2004). 



[34] S. Nojiri and S. D. Odintsov, arXiv:0804.3519 [hep-th]. 

[35] T. Koivisto, Class. Quant. Grav. 23, 4289 (2006). 

[36] I. Brevik and O. Gorbunova. larXiv:0806.13 99 [gr-qc]; to appear in Eur. 
Phys. J. C. 



12 



